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Abstract: The Magnetosplieric Multiscale (MMS) mission consists of four identically instrumented, 
spin- stabilized observatories, elliptically orbiting the Earth in a tetrahedron formation. For the 
operational success of the mission, on-board systems must be able to deliver high-precision orbital 
adjustment maneuvers. On MMS, this is accomplished using feedback from on-board star sensors in 
tandem with accelerometers whose measurements are dynamically corrected for errors associated 
with a spinning platform. In order to determine the required corrections to the measured accelera- 
tion, precise estimates of attitude, rate, and mass-properties are necessary. To this end, both an 
on-board and ground-based Multiplicative Extended Kalman Filter (MEKF) were formulated and 
implemented in order to estimate the dynamic and quasi- static properties of the spacecraft. 
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1 Introduction 

The Magnetospheric Multiscale (MMS) mission, launched on March 13, 2015, is the fourth mission 
of NASA’s Solar Terrestrial Probe program. MMS consists of four identically instrumented obser- 
vatories that function as a formation to provide the first definitive study of magnetic reconnection in 
space. 

Since it is frequently desirable to isolate electric and magnetic field sensors from stray effects caused 
by the spacecraft’s core-body, the suite of instruments on MMS includes six radial and two axial 
instrument -booms with deployed lengths ranging from 5-60 meters (see Fig. 1). The observatory is 
spin- stabilized about its positive z-axis with a nominal rate slightly above 3 rev/min (RPM). The 
spin is also used to maintain tension in the four radial wire-booms. 

Each observatory’s Attitude Control System (ACS) is responsible for orbital adjustments, attitude 
control, and spin adjustments. Its sensor and actuator compliment consists of Adcole’s spinning 
slit digital sun sensor, a four camera-head pASC Star Tracker System (STS) manufactured by 
the Technical University of Denmark, the micro-g resolution Acceleration Measurement System 
(AMS) produced by ZIN Technologies, and four axial AMPAC 1-lbf (4.4 N) and eight radial Aerojet 
4-lbf (17.8 N) mono-propellant hydrazine thrusters. Additional details on the ACS hardware and 
maneuvering system controller may be found in the references[l][2]. 

The sections that follow describe the formulation and performance of the Multiplicative Extended 
Kalman Filter (MEKF) that is used for on-orbit attitude and rate determination, and the addi- 
tional augmented filter-states and measurements necessary to perform an off-line mass property 
and thruster calibration. The paper concludes with a description of a calibration maneuver se- 
quence performed with the spacecraft in a partially-deployed state, and a comparison of the system 
identification results against pre-flight expected values. 
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Figure 1. MMS Observatory Fully-Deployed (top) and Stowed (bottom) 


2 Multiplicative Extended Kalman Filter (MEKF) 

The MEKF is an evolved version of the Extended Kalman Filer (EKF) that was originally applied 
to the Space Precision Attitude Reference System (a.k.a. SPARS) in 1969, and has been rigorously 
developed by Lefferts, Markley, and Shuster[3][4][5]. It has frequently been used for spacecraft 
attitude determination, in both real-time and off-line systems. In particular, the MEKF addresses 
a pair of difficulties associated with using a quaternion in an Extended Kalman Filter — namely, 
preserving the four-component quaternion’s unity norm constraint, while maintaining an unbiased 
estimate of the attitude. The fundamental idea of the MEKF is to use a quaternion as the “global” 
attitude, and use a three-component “local” representation for the attitude error (states) in the EKF. 
Whereas a more intuitive EKF approach might use all four elements of a error quaternion (Aq) as 
states, and then perform a full-state update additively by q true = Aq + q (followed by an ad hoc 
normalization scheme), the MEKF advocates a multiplicative full-state update 

qtme = 5q 0 q (1) 

where q trU e is the quaternion parameterization of the true rotation/transformation from an inertially- 
fixed reference frame to a body-frame attached to the spacecraft, q is the (unbiased) estimate of 
q trU e- and bq is an error quaternion that is parameterized by three states in the filter — the error 
vector, 86. Note that q tme , q, and <5q are all properly normalized unit quaternions. 
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Also introduced in Eq. (1) is the notation for a (left) quaternion product that denotes a 4 x 4 matrix 
formed from a four element column-vector, that expands to 
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with q 1 : 3 denoting the vector part of the quaternion and <74 the scalar element. This notation is 
analogous to the skew-symmetric matrix that is formed from a three-component base-vector, e.g. 
the angular rate vector uj, as 
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that is equivalent to the vector cross product in matrix algebra. Lastly, the symbol I„ is used through 
the text to represent the identity matrix of dimension n x n. 


What exactly then is the error vector SO that appears in the filter? Well, for small errors, a lst-order 
approximation is a legitimate abstraction — and the error states of a Kalman Filter are typically kept 
small. Conceptually, a perfectly adequate “mental picture” of 56 is a vector of small error angles 
(e.g. d x ,d y ,d z ). With that in mind, the lst-order relationship between the error vector and the error 
quaternion is 


<5q = qtrue(q ) 1 



(4) 


which is now not of unit norm. However, even in the Kalman Filter’s algorithm — one that is linear 
by design — a lst-order approximation is an insufficient description of the error vector given that 
the ultimate goal is to systematically address concerns surrounding quaternion normalization and 
biased-estimates. For that purpose, a more precise definition of the filter’s three-component error 
vector 56 is needed. Fortunately, many rigorously defined three-component attitude representations 
are available to choose from — the Euler rotation axis and Euler angle, the Gibbs vector (a.k.a. 
Rodrigues parameters ), modified Rodrigues parameters, Tait-Bryan angles (a.k.a. a 1-2-3 Euler 
sequence), etc. — each with their own benefits and pitfalls. Which particular parameterization 
is “best” is entirely a system design decision, and the references ([3], [4]) expound on the usage 
for each possible choice. In fact, the only places in the MEKF where the particular choice of 
parameterization (i.e. exact definition of 56 beyond lst-order) comes into play are the measurement 
residual computation and reset operation — both of which will be discussed shortly. 


For MMS, the Gibb’s vector (g) was selected to be the three-component parameterization of the 
attitude error. It is related to the error quaternion by 
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or its inverse 
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and to the error vector states by 


se = 2Sg 


(7) 


This error vector definition allows for a more precise re-statement of Eq. (4) as 
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The Gibbs vector parameterization is a mapping from four-dimension quaternion space to three- 
dimensional Euclidean space, and is infinity for 180-degree rotations (i.e. q 4 — 0). While not 
recommended as a global attitude representation, it has some nice properties when used in conjunc- 
tion with the MEKF. These benefits include: 


• free of singularities up to ± 180-degrees (error states will never be that large) 

• avoids accumulation of numerical errors in the full-state quaternion norm through an explicit 
normalization in the “reset” that is neither an ad hoc re-normalization operation, nor does it 
require transcendental function evaluations (see Eq. (8)) 

• largest possible 180° attitude errors map to infinity, so the representation is conceptually 
compatible with a Gaussian (or other random distribution) with infinitely long tails 

• observation model is insensitive to the sign ambiguity in the star camera’s output quaternion 

• although not unique to the Gibbs vector parameterization, the diagonals of the error covariance 
matrix ( P ) map directly to attitude error variance (a 2 ) 


To summarize, the attitude error-states (SO) used in the MMS MEKF are equal to twice the Gibbs 
error-vector (fig), and for “small” values are approximately the angle of rotation about each of the 
body-axes from the current full-state estimate (q) to the true attitude of the spacecraft (q tme ). 

2.1 Dynamical System Model (MMS ACS) 

Armed with a better understanding of the nuances of the MEKF, we may now focus on the specifics 
of the MMS application. The attitude control system (ACS) flight software (FSW) is limited 
to the computational capabilities of the flight processor (Motorola RH-CF5208 Coldfire), which 
lacks hardware acceleration of floating-point operations. As a result, the ACS FSW models the 
spacecraft’s angular rate (cu) dynamics using Euler’s rotational equation of a simple rigid-body 
(even after all the boom appendages are deployed) that is driven by autonomously commanded 
(u(?)) thruster torques (r(u)). 


r(u) = Icu +uj x luj 


(9) 
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where I is the second mass moment matrix about the spacecraft’s center of mass (r c ). The true 
attitude kinematics are in turn driven by the angular rate. 
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The commonly used gyro substitution model[4 ] is not an option for MMS due to the lack of an 
on-board gyroscope. Therefore, the angular rate (expressed in the body-frame) along with the 
attitude quaternion (from inertial to body-frame) make up the seven states of the dynamical system 
model that is the basis for the on-board MEKF. 


As previously discussed, the MEKF operates on only six states of attitude and rate errors — the 
number is reduced by one using the error vector SO attitude parameterization. Since the reset 
operation (discussed in section 2.3) moves the error state information into the full-state after 
each measurement is processed, and the full-states are propagated using (non-linear) Runge-Kutta 
integration, the MMS MEKF does not ever perform an explicit propagation of its error states 
(x). However, the error-state covariance (P) is propagated using linearized dynamics, so it is still 
necessary to determine the coefficients of the linearized (state-space) model. 


For the non-attitude error states (e.g. u >), the linearization process follows the standard EKF 
template[l 1]. The non-linear system dynamics (x = f(x,u,f)) are expanded about a reference 
trajectory (i.e. the current state estimate x(f)), using a first-order Taylor series expansion to obtain 
the error-state dynamics 
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In the case of the spacecraft angular rate error, we have 
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Recalling that the partial derivative of a scalar with respect to a vector is a vector, and that the 
partial derivative of a vector with respect to a vector is a matrix, we very much want to avoid taking 
the partial derivative of a matrix with respect to a vector (a hyper-matrix?)! With this in mind, a 
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matrix-equivalent to the vector cross product identity u x v = — v x u is used along with the chain 
rule to evaluate the partial derivative of Euler’s equation with respect to the angular rate — careful to 
always keep the vector being differentiated on the far-right of each compound term so that it drops 
out as the identity matrix. 
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which yields the following linearization for use in the angular rate error dynamics of Eq. (17), 
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(19) 


Unfortunately, the additive error definition (Sx = x — x) of Eq. (12) defeats the stated purpose of 
the MEKF (i.e. multiplicative update), and deriving the dynamics for the attitude-error requires 
traveling along a somewhat longer path. One route, taken by [4] and [5], begins by using the chain 
rule to obtain the time derivative of the attitude error definition of Eq. (1), 
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In Eq. (23) the quaternion identity p®q = q©p was used[4], which introduces here the “left” 
quaternion multiplication matrix notation, []®. The left-product is defined identically to the right- 
product (Eq, (2)) except for the sign on the cross-product term in the upper-left 3x3 block of the 
matrix. The expressions of Eqs. (22) - (24) for attitude error-quaternion are wonderfully exact. 
However there are two problems — the quaternion is not the final attitude-error parameterization 
that will be used in the filter, and it is non-linear in the error-states. In order to make the change of 
variables and linearize the error-state dynamics, the lst-order approximations <5qi : 3 « S6/2 and 
8q4 = 1 (see Eq. (4)) are used for the error-vector, and terms that are the product of two (small) 
error-states are dropped <8oj86 « 0). With these substitutions, Eq. (24) becomes 
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which is the desired final form of the error vector dynamics. The dynamics for the full and error 
states of the ACS flight software MEKF are summarized in Table 1. 



2.2 Measurement Update 

There are three sensors on the MMS observatories — a sun sensor (DSS), a star tracker (STS), and 
an accelerometer (AMS). Of the three, only the STS is used for on-board (real-time) attitude and 
rate determination. The DSS is excluded because its resolution of ±0.125° (450 arcsec) was shown 
to contribute little to the solution. The AMS acceleration measurements are neglected because of 
the modest capability of the flight processor (the filter would need to be augmented with bias states), 
and the additional power demands (AMS electronics plus thermal control heaters for bias stability). 
Nevertheless, the AMS is enabled for maneuvers, and the richness of information contained in 
its 1 kHz stream of \ig acceleration measurements makes it possible to perform off-line system 
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calibration using ground-telemetered data (the subject of the second-half of this paper). The STS 
and AMS measurements are staggered in time, so there are no update issues with simultaneous 
measurements. 


The specific form of the MEKF used on MMS is sometimes referred to as a Continuous-Discrete 
Extended Kalman Filter[6\, due to the discrete measurement updates from the on-board sensors 
combined with the continuous state-dynamics. The generalized nonlinear observation model is 

yk — h (qtrue (^) > k>true (tk) > * - -)jfc T W ( 26 ) 

where is a vector of Gaussian distributed random measurement errors with covariance matrix R k . 


The goal of the measurement update is to use feedback of the difference between the actual (noisy) 
measurement y/ t . and a prediction of what it “should be” (i.e. its expected value y k ) in order to adjust 
the estimate using an optimal feedback gain (K k ). In the parlance of a Kalman filter with discrete 
measurements, this is called the state update , and K k is the Kalman gain. This operation — which was 
described verbosely in the preceding paragraph — can be expressed more succinctly and elegantly 
with the following mathematical statement 


<5x^ = dx fc +K k {y k -y k } 


where the Kalman gain used is in the EKF standard form 
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Furthermore, in an EKF, the measurement model is allowed to be a nonlinear function of the states, 
and is historically segregated into two parts using a lst-order Taylor series approximation 
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where El k is known as the measurement sensitivity matrix — conceptually a local gradient in the 
nonlinear measurement with respect to a particular state. While this is pedestrian fare for a standard 
EKF, it has been made explicit here because it is another instance in the MEKF where the attitude 
error-state must be handled carefully. Specifically, for the additive states it is immaterial whether or 
not the sensitivity is with respect to the full or error state (i.e. the partial with respect to x produces 
the same result as the partial with respect to dx). Not so for the multiplicative attitude! It requires 
a more precise statement of the gradient. Since the linearized “slope” of the partial derivative 
multiplies the error states in Eq. (31), the measurement sensitivity matrix (H k ) is correctly defined 
with respect to error states. However, for the remainder of the paper the sensitivity partials of the 
additive states will still appear as dx (i.e. with respect to a full-state), solely because d (dx) is 
cumbersome notation. With that said, the definition of the measurement sensitivity matrix for the 
MMS on-board MEKF can be written as 


H k = 
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where it can be shown[4] by using the lst-order approximation of Eq. (4) produces the matrix 
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Finally, combining Eqs. (27) and (31) results in a state update for the kth measurement in the 
standard EKF form 
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where p k is referred to as the measurement residual. In the sections that follow, specific models are 
developed for the star tracker and accelerometer discrete measurement updates. 


2.2.1 Star Sensor Measurements 


The output of the MMS star tracker is simply a single quaternion per camera head unit (q c hu)- 
Because it is important that the measurement y k closely matches the estimated measurement y k in 
the residual calculation, the approximation of Eq. (33) is not used in favor of a measurement model 
that exactly matches our (Gibbs vector) parameterization of the attitude error. Specifically, the STS 
measurement model is 
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which possesses the previously mentioned benefit of insensitivity to sign ambiguity in the star 
tracker output. The associated measurement sensitivity matrix is then 
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and the star tracker measurement residual itself is 

(Pchu)fc = (ychu)^ h (q^ ) ^chu 


_l3 ®3x3 


q 


se k 

Sid, 


0 ((q C hu ,(Cq it 1 ) 1:3 rir ft l 

— 2 — r^r- — TT 2 , - ~ - — tt [JI 3 "3x3 J 


((qchu)®qjt ) 

= (^^chu)yt — 


(qf % 1 ) 4 


se k 

Sid 


k J 


•bxi 


(39) 


(40) 

(41) 

(42) 


Additionally, as will be explained shortly, the reset operation ensures 86 k is always zero. 
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For improved computational efficiency, the MMS ACS FSW combines the four simultaneous camera 
head unit observations into a single measurement to be processed by the MEKF[7]. This “effective” 
measurement is just the individual camera head residuals weighted by their associated measurement 
covariance matrix R" hu (transformed in to a common frame). This approach is tractable because all 
four measurements share a common sensitivity matrix H c h u . 

(EWiuC) (43) 

4 

= (Krff)t £ (^hu)i ' Wh„)t (44) 

n= 1 

On-board MMS, the measurement covariance matrix R^ — E {wvj} is calculated one of two ways. 
It is either a time- varying matrix that is constructed from the CHU’s internal algorithm’s reported 
image fit-residual (a single scalar unsigned byte), or else it is a constant (diagonal) matrix loaded into 
the FSW. The latter was originally based upon the DTU STS performance specification (a_ — 20, 
°boresite = 60 arcsec), and then later tuned to the statistics of the actual flight residuals. Of the two 
methods, the fixed R approached proved to be more accurate. 

The ensemble-average of the star camera head solutions essentially defines the spacecraft’s body- 
frame. The knowledge error with regard to the alignment of the heads (relative to their nominal 
design and verified by ground metrology) is handled by the separate off-line calibration system 
known as the MMS Attitude Ground System. The goal of that calibration is to minimize persistent 
biases from the measurement residuals by adjusting the alignment estimate of each CHU. Details of 
the calibration process can be found in [8]. 

2.2.2 Accelerometer Measurements 

The high-performance, tri-axial, acceleration measurement system (AMS) manufactured by ZIN 
technologies exists on the observatories in order to perform precise closed-loop orbital adjustments 
of the MMS formation) 1] [9]. However, the AMS’s high-rate sensor data is also available to perform 
off-line system identification. Without delving too deeply in the derivation provided in [1], we can 
import the final form of the sampled acceleration measurement model 

(yams)yt = hams ( f t h r 5 ^*C5 (45) 

a k = (r£ - r c ) +u£u£ (r d -r c ) - ^(2 • a^ x r c + f c )^ +b fc + (v ams ) t (46) 

r cd multi-body effects 

where m is the mass of the spacecraft, f t h r is the sum of the body-fixed thrust-force, r d is the location 
of the accelerometer in the body-frame, r c is the location of the combined center-of-mass of the 
entire observatory, b is a vector of the accelerometer’s intrinsic biases, and v am s is the measurement 
noise. For a discussion of the other potential measurement errors and their mitigation, see the 
reference[l]. 

A rigid-body assumption allows us to disregard center-of-mass motion (r c ,r c = 0), which results in 
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the accelerometer’s measurement sensitivity matrix with respect to the attitude and rate error states 
as 
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The expression for the measurement partial with respect to can be obtained from Eq. (46) using 
the chain rule and prudent cross product term-swapping. 
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The preceding result was obtained using Eq. (18) and the identity (a x b) x 


a x b x —b x a x . 


The covariance of the measurement noise can be derived from the AMS specification which asserts 
the root-mean-squared noise (a rm s ) should be less than 8 /igrms in the frequency range 0-10 Hz, and 
80 /igrms from 10-500 Hz. This implies that if a (single-sided) periodogram of acceleration-samples 
over a large time-interval was computed, its expected magnitude (A) over a given frequency range 
(e.g. A/h z = 10 Hz) should be less than 
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Figure 2 shows this is indeed the case for a sample taken from a non-maneuvering MMS-1. The 
AMS anti-alias filtering with break-point near 250 Hz can clearly be seen in the result. This implies 
a measurement covariance matrix with diagonal elements of R « 3.8 x 10 7 is reasonable, although 
the running sample-variance inside the AMS typically reads lower ( a ~ 16/ig) on most units. 



frequency [Hz] 

Figure 2. PSD of MMS-1 AMS 1 kHz Acceleration Data 
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2.2.3 Covariance Update 


The final element associated with the receipt of new measurement observations, is the update of 
the error-state covariance P. For the MMS MEKF, Jordan ’sform of the update equation is used to 
improve numerical stability 

P+ = (1 6 -K k H k )Pt(I 6 -K k H k ) T +K k R k Kj (51) 

and the resultant is forced symmetric using averaging (of round-off errors) 

p t=\(Pk+(K) T ) ( 52 ) 

To date, no stability problems have been experienced with the flight software MEKF that has 
operating continuously in single-precision floating-point for the past seven months. 

2.3 Reset 

Intrinsic to the MEKF is a process by which the information contained in the error states after a 
measurement fix f is transferred to the pre-update full state estimates x , while “simultaneously” 
resetting the error state to zero. For the attitude state — using a Gibbs vector error parameterization — 
the reset operation looks like 
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that effectively (and legitimately) enforces quaternion normalization. 80 is now also set to zero, 
and (as shown in [4]) stays zero throughout propagation. 


The reset of non-attitude states follow the standard linear addition update of an EKF. For example, 
the angular rate is updated by 


<JL> + — <JL> + 8<JL> + 
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and 8u> is set to zero. Because all of the reset/propagated error states are zero (SO ,8oj = 0), 
Eq. (34) collapses to simply 
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that — for what its worth — permits a so-called implicit [4J (i.e. combined) measurement/full-state 
update for the non-attitude state(s) 


+K COk [y k -h(q k ,u k )] 

effectively replacing Eqs. (34) and (57). 

2.4 Propagation 


(59) 


Certain estimated quantites of the MEKF must be propagated from one measurement to the next. 
For the MMS ACS that propagation interval is typically the period of its 4 Hz control cycle, 
which is broken in to two sub-intervals — the center-of-integration time stamp of the star tracker 
measurements, and the arrival of the AMS packet (in order to transform and utilize the incremental 
velocity for closed-loop maneuvering). As previously mentioned, the nonlinear models expressed 
as ordinary differential equation in section 2.1 governing the full-state dynamics are propagated 
using a Runge-Kutta integrator. Due to the reset operation, the error-states are zero and do not need 
to be propagated. However, the covariance of the error-state vector does. 

Pk +l =^kPk^ T k +Qk (60) 

where Qk — E { w/- wj } . This is done discretely using the linearized error-state dynamical models 
derived in section 2.1 and summarized in Table 1. The linear system dynamics of the estimate 

8x(t) — F(t)8x(t) + G(t)w(t) (61) 


are assumed constant over a small interval (At) so that the dynamics may be discretized as 


8x k+l =<Mx+ + w/, 


(62) 


and the derived state and noise-input transition matrices (<t>. Ej may be used to propagate the error 
covariance. The most exact route to discretized the error-state dynamics is by using the method of 
Van Loan [10] with the matrix exponential applied to the following construction 


A 


—F GW G J 

0 F J 


(63) 


and then extracting the desired sub-matrices needed for the covariance propagation from 


e AAf = 


B 

0 


*kQk 

* 2 l \ 


(64) 
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For the flight software, it was entirely impractical to consider using the exponential to construct the 
necessary matrices in real-time. Instead, the ACS FSW resorted to the lst-order approximation 

<S> k ^I 6 +F{t k )At (65) 


and also from [11] 

Mr 

Qk= d> fc (Ar-T)Wd>(Ar-T) T dT (66) 

Jo 

where W = diag{Wg, W m } is the diagonal matrix of spectral amplitudes (variances) for the “white” 
process-noise w (t) that drives the attitude and rate dynamics. Using the approximation for from 
Eq. (65) yields, 

lW 0 At+(W 0 ^-^W e )^ + W m ^- + W (0 F 0) ^- 

[ Wco^ + F^Wco^- W m At + (F m W a + W 0J Fj } )^+ F a W a fJ } 

(67) 

While the process noise is potentially very small, the MMS MEKF uses values around 1 0 6 to 
purposely “de-tune” the filter as protection against un-modeled dynamics, parameter errors, and 
maintain numerical stability. The matrix F a is the lower-right sub-matrix of F(t), and is given 
by Eq. (19). Additionally, the ACS FSW propagation of Eq. (60) also made use of sparse-matrix 
multiplication optimizations to help reduce some of the MEKF computational load. 

3 Augmented States for System Identification 

In addition to the on-board filter that has been described, a more elaborate variant was developed 
for ground calibration. This system-id MEKF adds twelve states for mass property estimation, and 
three additional states to determine the steady-state force due to each thruster. The motivation for 
the work will now be described. 

As a formation, the MMS mission must perform orbital maneuvers with an accuracy of roughly 1% 
(3c) or else potentially fall in to a “tail-spin” of continuous orbit corrections. In order to satisfy 
this accuracy requirement, an AMS was added to each of the four observatories, and algorithms 
developed for velocity control using real-time feedback from the accelerometers[l]. By going 
“closed-loop”, MMS successfully mitigated most of its parameter sensitivities (e.g. mass knowledge), 
while exposing itself to some others. Two system properties in particular were identified as the 
most potentially damaging while simultaneously the most difficult to estimate on the ground — the 
composite spacecraft center-of-mass (CM), and the steady-state thruster force. 

The CM is important because knowledge of it’s lever-arm to the accelerometer sensor-heads (r CY /) is 
a key component of the centripetal-compensation algorithm used in the ACS for inertial velocity 
tracking. This compensation attempts to mitigate the velocity-errors that occur due to the second 
and third right-side terms of Eq. (46) when it is numerically integrated. Additionally, knowledge of 
the CM location is used (on the ground) to determine thruster duty-cycles that effectively balances 
thrusters into usable pure-translations (i.e. “torque-less”) combinations. 
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The second-half of the “thrust-balancing” equation is the magnitude of the thrust itself. It was 
shown using Monte Carlo analysis, that steady-state thrust uncertainty of greater than ±3% (3 ( 7 ) 
could result in an imbalance that would unacceptably degrade maneuvering performance and/or 
exceed stress limits on the structure (e.g. the ADP axial booms). It is important to note that 
accurate thrust-knowledge is actually needed due to secondary effects on the system, because 
closed-loop control has already eliminated it’s primary effect. If that were not the case, a 3% error 
in thrust-knowledge would never permit a 1% maneuvering accuracy. Q.E.D. 

Lastly, the second mass moments of inertia were added to the filter (reluctantly) when it appeared 
that the estimation processes simply would not perform well without them. Inertia mismatch 
between the MEKF model and the plant-truth aliased itself as errors in many other state estimates, 
and prevented the filter from achieving the desired degree of accuracy. The specific form in which it 
will be introduced — using the fuel’s contribution as states, instead of using the entire observatory’s 
mass properties directly — was also a just a progression of the design based on trial and error that 
departs from [12]. No claim is being made at this time that it represents a fundamental characteristic 
of an effective estimator, or even that it improves convergence, stability, etc. With this addition, the 
total augmented error-state vector (assuming a single thruster) has become 

Sx— [SO Sul 5b 5rf 5If 5f w 5T C 5T x ] T (68) 


3.1 Accelerometer Biases 


In order to effectively estimate any of the system characteristics it is necessary to incorporate the 
accelerometer measurements. This can only be accomplished accurately if the accelerometers’ 
intrinsic (thermo-electric) biases b( true ) are also estimated by the filter. The dynamics model for the 
biases is simply 

b = ffc(x(f),w(f)) =0 + Wfc(f) (69) 


since we assume the three axial biases are effectively constant over the relatively short calibration 
period (< 2 hours), but could choose to set the process noise W/ ; to reflect a small bias drift — 
specified in the AMS to be less than 1 jig over 12-hours. The dynamic model is obviously identical 
for both the full-state dynamics (b) and the error states (5b), so the linearized error-state matrix 
expands to 
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(70) 


For a measurement y ams taken by the accelerometer (Eq. (45)), the associated measurement sensi- 
tivity with respect to the bias is 


dh ams 

db 


(71) 
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3.2 Fuel Mass Perturbations 


A careful bookkeeping and spin-balancing campaign[13] obtain what is believed to be very accurate 
knowledge of the stowed MMS observatories’ dry mass properties. However, tests that practiced the 
loading and unloading of the 41 1.6-kg of fuel from the stiff elastomeric diaphragms inside the four 
semi- spherical tanks showed that the fuel-mass could take on a variety of “shapes” which could in 
turn cause significant perturbations in the total system mass properties [14]. 

Nine states need to be added to the filter in order to estimate the effect of the fuel mass property 
uncertainty — three for the location of the combined (in the four tanks) center-of-mass of the fuel 
Tf in the body-frame, and six for the fuel-mass’s moments and products of inertia about it’s CM, 
If = [ifxx If If zz If If xz If yz ] T . Ground tests performed by Southwest Research Institute 
on a mock-up of a single MMS fuel tank, indicted that for a full tank, the mass participation in the 
fuel modes was relatively small («18%) and at a frequency of around 5 Hz[15]. Since this was not 
expected to significantly affect the system dynamics (and high-fidelity simulation of the truth model 
verified it doesn’t), the fuel state dynamics are also modeled as constants 

r f = f r (x(f),w(f)) = 0 + w r (f) (72) 

if = f,-(x(f),w(f)) = 0 + w,-(f) (73) 

The only non-zero partials in the linearized error-state dynamics are di' m /drf, and c)f 0) / r)If. Of 
the two, an analytic expression could only be found for the first. For the partial with respect to the 
“vectorized” fuel inertia, a numerical finite-difference solution was used to overcome the difficulties 
with taking partial derivative of a matrix with respect to a vector. While the finite-difference 
approach could also be brought to bear on the fuel center-of-mass partial derivative as well, the 
analytical solution is presented here because it is likely to be more computationally efficient, and 
because of the insights it can provide regarding the partitioning the system mass into dry and wet 
components. 


The composite observatory center-of-mass, r c , can be broken down into wet and dry components 


r c 


in dry 

m 


r dry + 


m fuel 

m 


(74) 


where the total mass of the spacecraft m = mf ue | + m dr y, and r c i ry is the location of the dry- structure’s 
center-of-mass in the body-frame. Utilizing a special case of the parallel axis theorem , the observa- 
tory’s moment of inertia about it’s center-of-mass, I, can be expressed as a sum of its inertia about 
the body-frame origin, J, and the “moment arm” from the origin to the center-of-mass (r c ), and then 
broken down into the fuel’s contribution and the dry-structure’s contribution. 


I = J + mr*r* (75) 

(jfuel T Jdry) mV c r c 

= (Ifuel -»hue|r f X r f X ) + (ldry-»Uryr d X ry r dry) +mT*r* (76) 

Note that in the present notation, If ue | is a 3 x 3 matrix where as If is a 6 x 1 vector. Proceeding 
with the derivation, Eq. (74) is combined with Eq. (76) to yield an expression for the spacecraft’s 
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moment of inertia as an explicit function of the new state rf 


I Ifuel T Idry T m dr 
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(77) 


Using an alternate form of the angular rate dynamics with the control-torque shown as a cross 
product of the thrust-force, and then taking its partial derivative with respect to the fuel CM yields 


1 lo+lo x 1 u = (r thr -r c ) x f thr 
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(78) 


(79) 


The terms needed to evaluate Eq. (79) can be determined by taking the partial derivatives of Eqs. (74) 
and (77). Specifically, they are 
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(81) 


The final required quantity, dld>/drf, is obtained by directly swapping uo with Co everywhere it 
appears in Eq. (81), and then estimating Co using the expected values of Eqs. (78) and (74) (i.e. by 
substituting in the current state estimates in place of truth-states in Euler’s equation). 


The acceleration measurement sensitivities associated with the new states are 

(82) 
(83) 

which may be evaluated using the results already calculated for Eq. (79), and the aforementioned 
finite-difference result for the fuel-inertial partials. 
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3.3 Thruster Dynamics 


The most difficult part of estimating the thruster output was due to “warm-up” effects. Based 
on ground-test data, it is expected that the MMS thrusters take nearly 15-seconds to reach true 
steady-state operation. It was unreasonable to expect that a single -thruster-at-a-time calibration 
(recall, the goal is the individual thruster modulation duty-cycles) could accommodate such a long 
period of “wasted” un-balanced thrust as part of its process. Instead, the data from the thruster 
delta-qualification program was mined in order to construct a dynamic model of the warm-up effect. 
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The result was a heuristic two-node thermal model, combined with a polynomial fit of the classical 
thrust-coefficient. The two thermal states along with a state representing the steady-state force itself, 
adds at least three states to the MEKF (but potentially three per thruster, depending on how the 
designer chooses to handle the individual thruster calibration in the filter software). 


The first component of the new model is a two node/state thermal model that approximates the 
heating and cooling of the thrust-chamber (where the catalyst-bed sits and two-stage hydrazine 
decomposition occurs). A two-node model was selected after the observation that there was no single 
cooling decay-constant capable of fitting the data “well”. The only elements of the models that are 
known from first-principles are: the flame-temperature of the hydrazine reaction (7f| ame = 1 875° F 
[1024 C]), and the cooling is either conductive/convective (linear in temperature) or radiative 
(quartic in temperature). A schematic of the model is shown in Fig. 3. 



heat of 
combustion 



chamber structure 

temperature temperature 

(state 1) (state 2) 

Figure 3. Thruster Thermal Node Model 


ambient 

temperature 


The full nonlinear (process) dynamic models for the two new thermal states are 

T c = f rc (x(f),u(t), w(f)) = u(t)-kf (7f| ame - T c ) + k x (T x — T c ) + w tc (t) (84) 

T x = f tt (x(f),u(f), w(f)) = k x (T c - T x ) + k (Too - T x )+k r (it ~ T?) +w tx (t ) (85) 

S V '' N v ' S V ' 

conductive radiative noise 

where u(t) is the control input (binary valve open-close status) for a particular thruster, and kf, k x , 
and k r are constant coefficients that were “hand-tuned” to fit the MMS qualification data (and later 
the flight data). The magnitude of thrust-force produced (/ t h r ) is found using a 3rd-order polynomial 
fit (/?o- Pi -Pi-.P'i) of the chamber temperature- state multiplied by the third state to be added to the 
filter — the steady-state force magnitude (f ss ) 

/thrO) = [p3T?(t)+p2T?(t)+piT c (t)+p 0 ] f ss (86) 

The thrust-vector — as used in Eq. (78) — is found by multiplying the magnitude by a unit-vector 
in the direction the thruster’s nozzle is pointing (d t hrX and regulating the output with the (binary) 
control input signal, u(t). In equation form, the preceding statement may be written as 

M0=/thr(0-dthr-w(0 (87) 

where the polynomials coefficients (/?„’ s) were again derived from the MMS qualification data — as 
was the somewhat arbitrary decisions to use a 3rd-order fit. 

Fastly, the true parameter of interest, the steady-state thrust-force f ss , is a function only of the 
hydrazine inlet-pressure in the propulsion system. Over the calibration period (e.g. a few short 
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thruster-pulses), this value is assumed to be constant. As has already been well established, the 
dynamic model of a constant, driven by process noise, is 

fss = /m(x(0>w(0) = 0 + Wf ss (t) (88) 

The effectiveness of this model in matching the qualification-test data is evidenced in Fig. 4 by 
how well dashed green line representing the simulated chamber temperature state (T c sim) tracks 
the blue line of measured chamber temperature data (T c meas), as well as the normalized force 
curves (black dashed-line vs. red measurement-dots) for a pulse-train of three consecutive valve 
cycles each continuously expelling 3.3-second of fuel on a 20-seconds interval (taken from a 
non-pulsed-width-modulated/class-A-pulse data set for the radial thrusters). 



Figure 4. Comparison of Simulated vs. Measured Temperatures and Forces 

The non-zero partial derivatives for the three new thruster-states, for use in the linearized error-state 
dynamics matrix F(t), are 
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The partial derivatives of the steady-state thrust force with respect to the thermal-states needed to 
complete the evaluation Eqs. (89) and (90) are 

^ = [3p 3 f c (t) 2 + 2p 2 f c (t) + Pl ] -f ss (95) 

rJ '( x,U 

=p 3 f c 3 (t)+p 2 f L 2 (t)+ Pl f c (t)+p 0 (96) 

Of™ x,u 

The non-zero measurement sensitivity partial derivatives with respect to the new thruster- states are, 
for the accelerometer measurement equation 

= -d thr — r^I -1 (r thr -? c ) x d t hr • u(t )-^f L (97) 

_ m J ujss x,u 

= -dthr-f^I _1 (r thr -f c ) x d thr ’ u(t) ■ (98) 

_ m J dl C x,U 

where the trailing terms — the thrust-force magnitude partials — are once again obtained by using 
Eqs. (95) and (96). 

4 Filter Performance 

Having defined the structure of both the ACS flight software MEKF and an augmented-state version 
for ground calibration, we now turn our attention to a battery of performance results for each. 
The two sources available by which to gauge performance are the MMS high-fidelity, nonlinear, 
time-domain simulation constructed in GSFC’s Freespace Simulation Environment[l6], and ground- 
processing of telemetered flight data from the observatories. 

4.1 Smoothing 

Before launching in to the performance result, one small aside is also worth mentioning. In order to 
access the true performance of the methods — and especially to obtain a definitive attitude and rate 
solution from the flight telemetry for which “truth” is not available — a Rauch-Tung-Stribel (RTS) 
optimal batch smoother was applied to the filtered results[17]. The RTS smoother falls in to the 
category of fixed-inten>al smoothers, and although it does not improve the final state-estimates from 
a filter operating on a finite sequence of data, it can significantly improve the estimates internal to 
the interval. 



The algorithm for the discrete-time RTS smoother uses stored values for x7, xf, fly.. P k , and P k 
from the forward-time MEKF to sweep backwards through the data and improve the estimates and 
error-covariance. The smoother initializes it’s recursion with the final (A-th) forward time estimates 
of the MEKF, x smN = x(; and / J sm v = P£, and proceeds (backwards in k) in the following manner [6] 


*sm k = *k +^sm k (x S m < . + 1 ~^k+l) 

Psm k — P£ ~ K smk (P k+l ~ Psm k+1 ) < 


(99) 

( 100 ) 

( 101 ) 
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The only downside to using the fixed-interval smoother is the burden of matrix storage, the fact 
that is it only possible as batch process, and that it does not improve the final estimate (for system 
identification). Nevertheless, many of the subsequent plots of the key system parameters-of-interest 
will also show an annotated curve for the smoothed estimates. 

4.2 Attitude and Rate Determination Performance 

4.2.1 Simulation 

A set of comparative plots shows the MEKF attitude and rate estimates surrounding a single 2.5- 
second pulse from MMS thruster #1. This unit is radially directed in the observatory’s minus-y 
direction (i.e. it applies a +y force), and is located approximately 0.5-meters below the CM. 
A “4-lb” thruster, at beginning-of-life pressure, applies 19.8 N of steady-state thrust. The mini- 
maneuver induces approximately 2.5° of nutation, and roughly the same amount of precession in 
the spacecraft’s total angular momentum. Figure 5 shows the effect of the single-thruster pulse on 
the observatory’s angular rate, as well as the 1 kHz acceleration at the sensor heads. 
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Figure 5. Simulated Single-Pulse Maneuver 
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The next series of plots show — for the same simulation set-up as above — the attitude and rate errors 
of the various estimators developed in the preceding sections. It is clear from Fig. 6 that the off-line 
MEKF outperforms the flight implementation, and unsurprisingly that the RTS smoother (applied 
to the augmented-state filter) outperforms both. 
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Figure 6. Comparison of Filtered Attitude and Rate Estimates of Simulated Maneuver 


There are a number of reasons why the higher-order filter produces a more accurate result. First, 
it is interesting to note that the augmented-state MEKF will — with identical tuning of the process 
and measurement noise statistics — very closely match the results of the six-state (flight software) 
MEKF as long as the acceleration measurements are excluded. Recall that the FSW MEKF was 
de-tuned as protection against un-modeled dynamics and numerical instability. With the acceleration 
measurements utilized by a filter with elevated attitude and rate process noise — while using an 
optimal measurement noise covariance — the relatively noisy AMS observations “drag the rate 
estimate about” too heavily. Either of two solutions will suffice to fix the issue — increase the noise 
statistics of the acceleration measurements, or (preferably) reduce the rate-dynamics process noise 
covariance. If either of these rectifications are applied, the results shown in Fig. 6 are achievable. 
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4.2.2 Flight Telemetry 


Figure 7 shows recorded telemetry from the ACS during a calibration maneuver performed on 
the observatory designated MMS1 that took place on April 1, 2015 03:12 (UTC). This maneuver 
(EA019) occurred while the observatories were in a semi-stowed state — only the ADP receiving 
element and the magnetometer booms were deployed. It was designed with the goal of exercising 
each individual thruster in matched pulses that induced — and then canceled — nutation. This was 
accomplished by spacing the 2.5-sec firings one -half of a nutation period apart (^15-sec). A couple 
of thrusters were fired four times in series to provide a more extended view of warm-up effects. 
Spin-rate change was essential for system-identification purposes — since it potentially separates the 
accelerometer-biases (static) from center-of-mass offset knowledge errors (a function of rate). Care 
was taken not to exceed the star tracker’s transverse or spin-rate limits, induce excessive pointing 
error (a thermal constraint), or significantly perturb the orbit (a collision risk). Nevertheless, it still 
made for a wild ride! 

These initial plots of the telemetry focus in on the first pulse in the maneuver. The telemetered 
data compares favorably with the simulation result shown in Fig. 5. The most glaring difference 


Telemetered Spin Rate Estimate AMS 1 kHz Structurally Filtered Meaurements 
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between the two is the acceleration plots on the right — the flight data is noticeably smoother. The 
reason behind this is that the raw-burst data from the AMS has undergone two operations in order 
to make it presentable — a temperature correction, and a structural (low-pass) filter. While the first is 
fairly benign, the structural filter eliminates almost all of the high-frequency noise characteristics of 
the AMS measurements. Unfortunately, this pre-processing step is needed because otherwise the 
effective rigid-body response is completely swamped by the excited structural modes. A high-order 
Finite-Impulse-Response filter with “brick- wall” frequency cut-off near 25 Hz was used because 
the linear phase-distortion (i.e. fixed sample lag) that it introduces is trivial to undo. 


A good metric for the performance of the flight MEKF is the measurement residuals for each of 
the star tracker’s camera head units (CHU) — shown in Fig. 8. They indicate that the ACS MEKF 
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Figure 8. ACS MEKF CHU Measurement Residuals During EA019 Maneuver 

is well-behaved — despite the aggressive maneuvering — and that the measurements are reasonably 
unbiased due to successful star sensor alignment calibration campaign[8]. Unfortunately, the 
statistics also show that the z-axis (bore-site) residual seems to be larger than expected. At roughly 
170 arcsec (la), the STS solution error is at nearly three times the spec-level of 60 arcsec (la). It 
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has been independently verified[8] that these elevated noise-levels about the bore-site persist during 
coasting — confirming that the phenomena is not an artifact of the ACS MEKF design, and most 
likely are just intrinsic to the tracker measurements themselves. 

4.3 Mass Property and Thrust Estimation Performance 

The effectiveness of the augmented-state MEKF at estimating the MMS mass properties and steady- 
state thrust magnitude for will now be demonstrated. As has already been asserted, the system-id 
filter performs nearly identical to the FSW MEKF when fed only star sensor measurements, and 
de-tuned to an equivalent degree. Despite have a significantly larger number of states, the system-id 
MEKF remains stable and converges in a reasonable vicinity of the true plant parameters. It also is 
well-behaved when given acceleration measurements at any interval between 0-1000 Hz. Despite the 
author’s concerns during the state-augmentation process, there appeared to be no serious information 
dilution pitfalls. 


One caveat regarding the augmented filter’s performance at estimating thrust is a sensitivity to 
correct timing of the control input, u(t). If the valve opening and closing times are not well known, 
or there is significant transport lag, the filter will tend to underestimate the thrust-level and/or over 
estimate the inertia. This touches upon one of the two coupled (aliasing) problems for MMS system 
identification. The thrust/inertia ambiguity is the first, and the center-of-mass/bias equivalency is 
another. Only through sufficient variance of excitation (i.e. spin-rate, nutation, thrust) can these 
quantities be resolved. Examination of the complete linearized measurement matrix H(t) 
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might offer some additional insight as to the observability of the augmented error-states (and which 
states need to be excited), but that analytical effort has not been pursued. Instead, evidence of the 
filter’s efficacy will be shown by way of simulation and Monte Carlo analysis. 


4.3.1 Simulation 

A computer simulation was performed that uti- 
lized a pair of thruster-pulses to perturb the 
spacecraft in a manner similar to the single- 
pulse of section 4.2.1. As part of the system- id 
post-processing, the errors listed in Table 2 in 
the initial estimates of the augmented (quasi- 
constant) states were introduced to demonstrate 
performance. 


Table 2. Summary of Initial Estimate Errors 


Parameter 

Error 

accelerometer biases 

+20 jig 

fuel center-of-mass 

+50 mm 

fuel moments of inertia 

+10 kg-m 2 

steady-state thrust magnitude 

+5% 


The process noise variances (i.e. the diagonal elements of the W matrix) were selected to be 1 0 1 6 
for the bias (b) and fuel states (rf. If). The variances of the thruster-related states (f s s, T c , T x ) were 
set larger (at 10 9 ) to account for model uncertainty. The results of the estimation process with 4 Hz 
star sensor measurements, and 100 Hz (decimated) acceleration samples are shown in Figs. 9-12. 
The high-fidelity simulation of the spacecraft truth-model with which the system-id MEKF is being 
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tested includes: colored sensor noise, time-varying mass (fuel depletion), fuel-slosh dynamics, and 
the thermal-dynamics of the thruster warm-up model — including small random actuation delays. 



time [sec] 

Figure 9. Observatory Center-of-Mass Estimation (Two-pulse sim) 



time [sec] 

Figure 10. AMS Bias Estimation (Two-pulse sim) 

Note that the results shown in Fig. 9 are for the total observatory center-of-mass — not just the 
fuel-mass that are the filter-states. Since the fuel-estimate has to potentially compensate for a 
number of parametric errors in the system, the total CM (and inertias) tend to be the more physically 
accurate result. 
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The bias estimates (Fig. 10) are also typical for MMS. The z-axis estimate tends to be the most 
accurate (which is actually the best possible result for the MMS application!!]), and the y-axis 
bias-estimate tends to be the least — most likely this is an artifact of the AMS placement relative to 
the spin-axis and center-of-mass. Nevertheless, bias estimation is not the primary goal of this filter, 
and the “tens of micro-g” error-regime of the estimate is sufficient for good system-id. 




0 10 20 30 40 50 60 70 

time [sec] 






Figure 11. Observatory Moments and Products of Inertia Estimation (Two-pulse sim) 

Figure 1 1 shows the evolution of inertia estimation process. This test-case shows that the 7 XZ and 
products of inertia are especially easy to identify on a spinner — a fact utilized on MMS to routinely 
perform a principal-axis (re)alignment as part of its pre-maneuver preparations!!]. 
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Figure 12. Thruster Steady-State Force Estimation (Two-pulse sim) 


4.3.2 Flight Results 

The angular rate estimates from the MMS1 ACS FSW for the full EA019 calibration maneuver are 
shown in Fig. 13. The stored real-time measurement data from the maneuver was processed by the 
off-line system-id MEKF in order to improve knowledge of the spacecraft’s operational parameters. 
The new calibrated mass properties, and derived thruster duty-cycles were then uploaded into the 
ACS flight software. The process was repeated for all four observatories, which have been operating 


ACS FSW AD (Principal Axis) Angular Rate Estimates 




Figure 13. MMS1 ACS Estimated Angular Rate Telemetry (Full EA019 Maneuver) 
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successfully [9] using the updated coefficients since April 2015. 


An example of the differences between the MMS 1 pre-flight and post-calibration parameter estimates 
are given in Table 3. The excellent agreement attests to the effectiveness of the augmented- state 
filter, and the careful accounting of the spacecraft assembly process [13]. 


Table 3. MMS1 EA019 Calibration Results 


Parameter 

Units 

Pre-launch 

Post-Cal 

Difference 

CM-x 

mm 

2.36 

2.86 

0.494 


CM-y 

mm 

5.03 

4.96 

-0.071 


CM-z 

mm 

607.68 

606.43 

-1 .244 


Ixx 

kg-m 2 

977.15 

975.33 

-1.82 

(-0.2%) 

lyy 

kg-m 2 

921 .24 

923.38 

2.14 

( 0.2%) 

Izz 

kg-m 2 

1594.99 

1594.84 

-0.15 

(-0.0%) 

Ixy 

kg-m 2 

-110.96 

-106.48 

4.47 

(-4.0%) 

Ixz 

kg-m 2 

-0.11 

-0.09 

0.02 


lyz 

kg-m 2 

-0.25 

-0.20 

0.04 


Thrust-01 

N 

18.423 

18.451 


1 .0799% 

Thrust-02 

N 

18.286 

18.334 


1.0731% 

Thrust-03 

N 

18.385 

18.371 


1 .0753% 

Thrust-04 

N 

18.329 

18.321 


1 .0723% 

Thrust-05 

N 

18.664 

18.635 


1 .0907% 

Thrust-06 

N 

18.727 

18.712 


1 .0952% 

Thrust-07 

N 

18.551 

18.515 


1 .0837% 

Thrust-08 

N 

18.425 

18.385 


1.0761% 

Thrust-09 

N 

3.901 

3.909 


0.9153% 

Thrust-10 

N 

3.926 

3.938 


0.9220% 

Thrust-11 

N 

3.798 

3.804 


0.8907% 

Thrust-12 

N 

3.868 

3.879 


0.9081% 


5 Multi-body Dynamics 

All of the filter development and flight results presented thus far have been based upon a rigid-body 
model of a spacecraft. For MMS, this was only a credible approximation during the commissioning 
period before the 60-meter long SDP wire -booms were deployed. Fortunately, the system identifica- 
tion maneuvers used to calibrate the mass properties and thrust were performed with the spacecraft 
in a semi-rigid configuration. 

In fact, by using a de-tuned MEKF and a set of “effective” mass properties (extrapolated from our 
semi-stowed calibration results), the ACS flight software is able to continue to use its six-state 
MEKF otherwise unaltered. However, details of the multi-body dynamics extend beyond the scope 
of the present text. The topic of this estimator’s performance with significant model-mismatch may 
instead become the focus of a future published work. 
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6 Summary 


In this text, some of the fundamental principals regarding the use of a Multiplicative Extended 
Kalman Filter were re-introduced in the context of spacecraft attitude and rate determination for the 
MMS mission. The use of gyroscopes and the gyro substitution simplification in Extended Kalman 
Filters has become so ubiquitous, that this application stands out for its absence. The flight version 
of MMS ’s six-state attitude and rate filter — operating solely off star camera measurements — was 
then documented in full, with the hope that the (admittedly pedantic) exposition might help clarify 
some of the conceptual nuances surrounding the MEKF. At the very least, this paper hopes to call 
attention to the excellent reference texts and papers that have recently become available, while 
evangelizing on the conceptual “tidiness” of the MEKF formulation. 

Additionally, the requirement for precise orbital maneuvering of the MMS formation belies the 
notion of a “simple spinner”, and drives the need for precise knowledge of the observatories’ mass 
properties and thrust-capabilities. A new estimator for system identification was presented that 
leveraged the flight MEKF framework — adding only high-rate accelerometer measurements and 
the desired quasi-constants as filter- states. With MEMS accelerometers becoming a fairly common 
component of a new generation of “CubeSats”, the hope is that this off-line process proves a useful 
weapon in slaying the dragon of system identification that can torture a mission designer. 

Finally, examples of the practical utility of the two filters was demonstrated using computer 
generated, and flight telemetered data. MMS’s unqualified success with the flight filter — even using 
a relatively low-powered processor (by today’s standards) — should further encourage the adoption 
of the MEKF as an aerospace staple. 
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